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Using samples of 2.25 x 10® J/^p events and 1.06 x 10* tp' events collected with the BES III 
detector, we study the /o(980) — >■ ao(980) and ao(980) /o(980) transitions in the processes J/ip 
(/)/o(980) (?!-a|^(980) and xa 7r"a|J(980) 7r"/o(980), respectively. Evidence for /o(980) 
Oo(980) is found with a significance of 3.4 a, while in the case of ao(980) /o(980) transition, the 
significance is 1.9 a. Measurements and upper limits of both branching ratios and mixing intensities 
are determined. 

PACS numbers: 25.75.Gz, 14.40.Df, 12.38.Mh 



I. INTRODUCTION 



The nature of the scalar mesons ao(980) and /o(980) has been a hot topic in light hadron physics for many years. 
These two states, with similar masses but different decay modes, are difficult to accommodate in the constituent 
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quark-antiquark scenario. Tremendous efforts in both experiment and theory have been made in order to understand 
them. Suggestions for their being exotic candidates, such as tetra-quark states, hybrids or KK molecules, can be 
found in the literature [1-6]. 

The possibility of mixing between aQ(980) and /o(980) was suggested long ago, and its measurement will shed light 
on the nature of these two resonances [7-14]. In particular, the leading contribution to the isospin- violating mixing 
transition amplitudes for /o(980) — ^ ao(980) and ag(980) — s- /o(980), is shown to be dominated by the difference of 
the unitarity cut which arises from the mass difference between the charged and neutral KK pairs. As a consequence, 
a narrow peak of about 8 MeV/c^ is predicted between the charged and neutral KK thresholds [8, 13, 14]. 

The mixing amplitudes strongly depend on the couplings of ag(980) and /o(980) to KK, and to probe the properties 

of these two scalar states, precise measurements of the mixing transitions are very important. Two kinds of mixing 

intensities, i.e. ^fa and ^af for the /o(980) — > aQ(980) and aQ(980) /o(980) transitions, respectively, can be defined 

and are accessible to measurement in charmonium decays [13, 14]: 

_ Br(J/V' 0/o(98O) ^ 0a;](98O) ^ 07?7rO) 
" Br{J/i} -> 0/o(98O) -> 4)7171) ' 

and 

_ Brixci 7r"ag(980) ^ 7r"/o(980) ^ 7r"7r+7r-) 
^ Br{xci -> 7r"aS](980) ^ ttOttO??) ' 

Using samples of 2.25 x 10** J/ip events [15] and 1.06 x 10^ V' events [16] collected with the BES III detector in 2009, 
we perform direct measurements of the aQ(980) — /o(980) mixing intensities via the processes J/t/i — > 0/o(98O) 
</>ag(980) ^ 07?7r" and xa ^ 7r"aSj(980) ^ 7r"/o(980) ^ 7:^71+71'. 

II. DETECTOR AND MONTE CARLO SIMULATION 

BEPC II is a double-ring e+e~ collider designed to provide a peak luminosity of 10^'^ cm~^s~^ at a beam current 
of 0.93 A. The BES III detector has a geometrical acceptance of 93% of 47r and has four main components: (1) A 
small-cell, helium-based (40% He, 60% CsHg) main drift chamber (MDC) with 43 layers providing an average single- 
hit resolution of 135 /xm, charged-particle momentum resolution in a 1 T magnetic field of 0.5% at 1 GeV/c, and 
the dE/dx resolution that is better than 6%. (2) An electromagnetic calorimeter (EMC) consisting of 6240 CsI(Tl) 
crystals in a cylindrical structure (barrel) and two endcaps. The energy resolution at 1.0 GeV/c is 2.5% (5%) in the 
barrel (endcaps), and the position resolution is 6 mm (9 mm) in the barrel (endcaps). (3) A time-of-flight system 
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(TOF) constructed of 5-cm-thick plastic scintillators, with 176 detectors of 2.4 m length in two layers in the barrel 
and 96 fan-shaped detectors in the endcaps. The barrel (endcap) time resolution of 80 ps (110 ps) provides 2(t Kj-K 
separation for momenta up to ^ 1.0 GeV/c. (4) The muon system (MUG) consists of 1000 m^ of Resistive Plate 
Chambers (RPCs) in nine barrel and eight endcap layers and provides 2 cm position resolution. 

The efficiency for J/ij] ^ 0/o(98O) — > 0aQ(98O) -> 0?77r° is estimated using a Monte Carlo (MC) simulation of 
J/'i/; (j)S, where S is the mixing signal represented by a narrow scalar Breit-Wigner uniformly decaying to rji:^ . The 
mass of the mixing signal is set to be 991.3 MeV/c^ (the center of {mj^+ +mx- ) and (m^'o +m|fo) [17]) and the width 
of the mixing signal is set to be 8 MeV/c^. The efficiency for V' IXa — > 77r°ao(980) — > 77r°/o(98G) — > 77r°7r+7r^ 
is estimated using a Monte Carlo simulation of ^Ij' — > 7Xci — > IT^'^ where S is the mixing signal with parameters as 
above, and decays into t:^t:~ isotropically. 

III. EVENT SELECTION 

Tracks of charged particles in BES III are reconstructed from MDC hits. We select tracks within ±20 cm of the 
interaction point in the beam direction and within 2 cm in the plane perpendicular to the beam. The TOF and 
dE/dx information are combined to form particle identification (PID) confidence levels for the tt. A', p hypotheses; 
each track is assigned to the particle type that corresponds to the hypothesis with the highest confidence level. 

Photon candidates are reconstructed by clustering EMC crystal energies. Efficiency and energy resolution are 
improved by including energy deposits in nearby TOF counters. The minimum energy is 25 MeV for barrel showers 
(I cos 6*1 < 0.80) and 50 MeV for endcap showers (0.86 < | cos 6*1 < 0.92). To exclude showers from charged particles, the 
angle between the nearest charged track and the shower must be greater than 10°. EMC cluster timing requirements 
suppress electronic noise and energy deposits unrelated to the event. 

The 7r° — >■ 77 and 77 — >^ 77 candidates are formed from pairs of photon candidates that are kinematically fit to the 
known resonance masses, and the from the kinematic fit with one degree of freedom is required to be less than 25. 
The decay angle of a photon is the polar angle measured in the 77 or tt*^ rest frame with respect to the 77 or tt*^ direction 
in the J/?/; or ?/;' rest frame. Real 77 and tt*^ mesons decay isotropically, and their angular distributions are fiat. 
However, the 77 and tt*^ candidates which originate from a wrong photon combination do not have a fiat distribution 
in this variable. To remove wrong photon combinations, the decay angle is required to satisfy | co&9decay \ < 0.95. 
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IV. MEASUREMENT OF J/tp (^/o(980) (^ag(980) (^7?7r° 

Events with two oppositely charged tracks identified as kaons and at least one distinct tt^ candidate and ry candidate 
are selected. A six-constraint kinematic fit (6C) is performed to the J/ip — ?► K^K^rjir^ hypothesis (constraints are the 
4- momentum of the J /tp and the tt" and rj masses), and x\c '^'^ required. If there is more than one combination, 
the combination with the smallest Xgc retained. The events are fitted to J/ip ^ K^K^rj-q and J /ijj — J^+lv ^7r°7r°, 
and the probabilities are required to be less than that from the kinematic fit to the signal channel J /ip K^rj-K^ . 

The K'^ K~ invariant mass distribution of selected events is shown in Fig. 1 (a), where a clear signal can be seen. 
The solid and dashed arrows show the signal and sideband regions, respectively. 

Figure 1 (b) shows the ryTr*^ invariant mass distribution recoiling against the 4> signal — 1.02| < 

0.015 GeV/c^). A narrow structure appears around 980 MeV/c^. The shaded histogram shows the rjn^ invari- 
ant mass of events recoiling against the sideband (1.065 GeV/c^ < Mj^+x- < 1.095 GeV/c^), and no sign of a peak 
near 980 MeV/c^ is evident. 

Exclusive MC samples of J/ip decays which have similar final states are generated to check whether a peak near 
980 MeV/c^ can be produced in the r/Tr" mass spectrum. The main backgrounds come from: (1) J /if.' — >■ ^tt^tt" 
via /o(980), /2(1270) or a phase space process; and (2) J/i' K*K-q + c.c, J/i^ K^il^'iQ)^!^ + c.c, J/ip 
7/2(1950) ^K*K* . The TyTr" invariant mass distribution from all these possible background channels is smooth, 
shown as the shaded region in Fig. 1 (c), and will not affect the determination of the number of mixing events. 

A MC sample of 2.0 x 10® inclusive J/il} decay events is used to investigate other possible backgrounds too. The 
shaded histogram in Fig. 1 (d) shows the "qir^ invariant mass distributions of events selected from the inclusive MC 
sample. In the 980 MeV/c^ region, there is no peaking background. 

An underlying process is from J/i/' (/)aQ(980) via a 7* or K*K loops [13]. However, it will produce a much 
broader distribution (50 — 100 MeV/c^) in the 7771° mass spectrum than /o(980) — > aQ(980) mixing [13]. We estimate 
its contribution in the fit to the mass spectrum. 

A simultaneous unbinned maximum likelihood fit to the rji:^ mass spectrums recoiling against the (j) and the 4> 
sideband is performed. In the signal region, the probability density function (PDF) is composed of the mixing signal, 
represented by the shape extracted from MC simulation of a narrow Breit-Wigner, the ao(980) contribution from 7* 
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FIG. 1: (a) The invariant mass spectrum of in J/ip — i- K'^ K" ■qiv'^ . The sohd arrows show the (j) mass window. The 

dashed arrows show the cj> sideband region used to estimate backgrounds. The ryvr" invariant mass of selected events is shown 
in (b), (c) and (d). The dots with error bars show the mass spectrum of M^^o recoiling against the (j). The shaded histogram 
is: (b) recoiling against the (j> sideband; (c) from exclusive MC; (d) from inclusive MC. 

or K*K loops represented by a Flatte formula/ and a 2nd order polynomial for the backgrounds. In the sideband 
region, the PDF is a 2nd order polynomial only. The mass of the mixing signal is set to 991.3 MeV/c^ (the center of 
(rnx+ +1TIK-) E^nd (rrixo +mj^o) [17]) and the width of the mixing signal is set to 8 McV/c^. The shape parameters 



^ The Flattc formula is taken from Rcf. [18]. The values of coupling constants and mass as used here are taken from the Crystal Barrel 
experiment results in Rcf. [19] and Rcf. [13]. 



of the background polynomials in the signal region and the sideband region are constrained to vary simultaneously 
in the fit. The normalization of each component is allowed to float. Figures 2 (a) and (b) show the results of the 
simultaneous fit. 

The fit yields N{fo a") = 25.8 ± 8.6 events for the mixing signal, and iV(ag (980)) = 13.6 ± 24.8 events for the 
ag(980) contribution from 7* or K* K loops. Comparing with the fit result without the mixing signal, the change in 
—InL with A(d.o./.) = 1 is 5.90, corresponding to a statistical significance of 3.4cr. Comparing with the fit result 
without the ag(980) contribution from 7* or K*K loops, the change in —InL with A((i.o./.) = 1 is 0.14, corresponding 
to a statistical significance of O.Scr. Using the Bayesian method, the upper limit for the number of mixing events is 
37.5, and the upper limit of the number of ao(980) from 7* or K*K loops is 50.8 events at the 90% confidence level 
(C.L.). The results are listed in Table I. 

The branching ratio of the mixing signal J/i/; — > (/)/o(980) — > 0ao(98O) 07^7r° is calculated as : 

Br{J/^ -> 0/o(98O) ^ K(980) ^ 0r?^°) = „ u ^ k+k-C r] ^ ^ n ( 0^ V 

£ fa ■ Nj/^ ■ Br{(p ^ K + A ) ■ Br{ri ^ 77) • Br (tt" 77) 

where iVj/^ is the total number of J/ip events and Sfa — (18.5 ± 0.2)% is the efhciency for the mixing signal 
J/ip (^/o(980) — > 0aj](98O) — > (f»-/Tr". The branching ratio is then determined to be (3.3 ± 1.1) x 10~^, where the 
error is statistical only. 

The total branching ratio of J/iJj -> 0aQ(98O) — >■ ^jyTr" is calculated as : 

BriJ/^P ^ 0,.«) = Ni.h^aO)/ef. + NiaO)/e. 

^ ^ ' ' Nj/^ ■ Br{(t> K+K-) ■ Br{T] 77) ■ Br(7r" ^ 77) ' 

where = (18.3 ± 0.2)% is the efficiency for the underlying process J/ip — > 0aQ(98O) — > (prjn^ . The branching ratio 
is then determined to be (5.0 ± 2.7) x 10^^, where the error is statistical only. 

If we fit the rji:^ invariant mass spectrum only with the mixing signal plus a 2nd order polynomial background, the 
fit yields 28.6 ± 7.0 events for the mixing signal. Comparing with the fit result with only the 2nd order polynomial, 
the change in —InL with A{d.o.f.) = 1 is 16.65, corresponding to a statistical significance of 5.8(t. The upper limit 
at the 90% C.L. is 39.1 events. 

If we assume there is no mixing and fit the TyTr" invariant mass spectrum only with the aQ(980) contribution from 
7* or K*K loops plus a 2nd order polynomial, the fit yields 75.8 ± 17.3 events for the ao(980) contribution from 7* or 
K* K loops. Comparing with the fit result with only the 2nd order polynomial, the change in —InL with A{d.o.f.) = 1 
is 10.89, corresponding to a statistical significance of 4.7cr. The upper limit at the 90% C.L. is 99.1 events. The results 
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TABLE I: Results of fits to the mass spectrum of ryvr' . The fit is described in the text, and the yield error is statistical 
only. N(/o — > ag) is the number of mixing events /o(980) — >■ ao(980)- S(/o — >■ flo) is the significance of the mixing signal 
/o(980) ag(980). N(a8(980)) is the number of aEi(980) events from 7* or K'K loops. S(ag(980)) is the significance of the 
ao(980) contribution from 7* or K* K loops. 
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FIG. 2: Results of the simultaneous fit of the t^tt'^ mass spectra: (a) recoiling against the 0, (b) recoiling against the 4> sideband. 
The solid curve in the result of the fit described in the text. The dotted curve is the mixing signal, and the dash-dotted curve 
is the ao(980) contribution from 7* or K* K loops. The dashed curves in (a) and (b) denote the background polynomial. 

are listed in Table I. The total branching ratio of J/i/j (/)aQ(980) — )• ^jyTr" is calculated to be (9.7 ± 2.2) x 10~^, 
where the error is statistical only. 



V. MEASUREMENT OF V' ^ 7Xci ^ 77r"ag(980) 77rVo(980) 77r"7r+7r- 



Events with two oppositely charged tracks identified as pions and at least three photons, which form at least one 
distinct 7r° candidate, are selected. A 5C kinematic fit is performed to the ip' 77r''7r+7r^ hypothesis (constraints 
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are the 4-momcntum of ip' and the tt^ mass) and xic ^ required. If there is more than one combination, the 
combination with the smallest xic retained. The events are also fitted to t/j' tt^tt+tt^ and -0' — > 7r°7r°7r+7r~, 
and the probabilities are required to be less than that from the kinematic fit to the signal channel ip' — > 77r°7r+7r^. 
To remove backgrounds with a J/tA decaying to leptons, the angle between the two charged tracks is required to be 
less than 160°. We further require the ratio of energy deposited by each charged track in the EMC to its momentum 
measured in the MDC to be less than 0.85. To remove the backgrounds which have jjJ/^ final states, the mass 
recoiling from any photon pair must not be in the J/ip mass window (|A^7ecoiimg ~ 3.097 GeV/c^| > 0.06 GeV/c^). 
The invariant mass distribution of tt^'tt+tt^ of the selected events is shown in Fig. 3 (a). 

The TT+TT" invariant mass distribution in the Xci mass window (3.49 GeV/c^ < M^o„+„- < 3.54 GeV/c^) is shown 
in Fig. 3 (b). A narrow structure around 980 MeV/c^ is evident. The shaded histogram shows the tt+tt" invariant mass 
of events in the xa sideband (3.39 GeV/c^ < M^o^+^- < 3.45 GeV/c^ and 3.54 GeV/c^ < A/^o^+^- < 3.59 GeV/c^). 

Exclusive MC samples of ip' decays which have similar final states are generated to check whether a peak near 
980 MeV/c^ can be produced in the tt+tt^ mass spectrum. The main possible backgrounds come from: ip' — > 'yXci 
^^J/ip — !► 777r"*'7r~7r*'; ip' 7r+7r~7r'^; ip' — > 7Xci ^ 7a^(980)7r^ — jrjir'^TT^T:^ and ip' — > rj'K^TT^'K^. The tt+tt^ 
invariant mass distributions from all these background channels is shown as the shaded histogram in Fig. 3 (c), and 
there is no peak around 980 MeV/c^. 

A MC sample of 1.0 x 10* inclusive ip' decay events is used to investigate other possible backgrounds. The shaded 
area in Fig. 3 (d) shows the 'k^tt" invariant mass distribution of events selected from the inclusive MC sample. In 
the 980 MeV/c^ region, there is no peaking background. The /o(980) from other ip' — ^ 7Xci -> 77r"/o(980) processes 
is much broader than the aQ(980) — > /o(980) mixing signal [14] and is estimated from the fit. 

A simultaneous fit is performed to the tt+tt^ invariant mass spectrum in the Xci mass window and the Xci sideband 
in a similar manner as in Sect. IV. The /o(980) contribution from other processes is represented by a Flattc formula^. 
Figures 4 (a) and (b) show the results of the simultaneous fit. 

The fit yields iV(a[J /o) = 6.4 ±3.2 events for the mixing signal and iV(/o(980)) = 0.0 ±8.6 events for the /o(980) 
contribution from other processes. Comparing with the fit result without the mixing signal, the change in —InL with 
A{d.o.f.) = 1 is 1.79, corresponding to a statistical significance of 1.9cr. Comparing with the fit result without the 



^ The Flatte formula is quoted from Rcf. [18]. The values of eoupling constants and mass used here are quoted from BES II experiment 
results [14, 20]. 
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FIG. 3: (a) The invariant mass spectrum of tt^tt^tt^ in ip' — > 77r+7r^7r''. The solid arrows show xa mass window, and the 
dashed arrows show the xa sideband region used to estimate background. The invariant mass of 7r"^7r~ of selected events is 
shown in (b), (c) and (d). The dots with error bars show the mass spectrum of M^+^- in the Xci mass window. The shaded 
histogram is: (b) Xci sideband; (c) from exclusive MC; and (d) from inclusive MC. 

/o(980) contribution from other processes, the change in —InL with A{d.o.f.) = 1 is less than 0.01, corresponding to 
a statistical significance of less than O.lcr. Using the Bayesian method, the upper limit for the number of the mixing 
events is 11.9, and the upper limit for the number of the /o(980) events from other processes is 16.7 events at the 
90% C.L. The results are listed in Table II. 
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The branching ratio of the mixing signal ip — > 7Xci — > 77r°aQ(980) — > 77r°/o(980) jtt^tt~^tt is calculated as: 



Br{ij' 7Xci ^ 77r°a[](980) ^ 77r°/o(980) ^ 77r''7r+7r~) 



iVK ^ /o) 



• iV^, ■ Br(7rO ^ 77) ' 

where N^f,/ is the total number of tp' events and Saf = (22.3 ±0.2)% is the efficiency for the mixing signal tp' 7Xci 
77r°a[](980) 77r°/o(980) — )• 77r°7r+7r^. The branching ratio is then determined to be (2.7 ± 1.4) x 10^'', where the 
error is statistical only. 

The total branching ratio of t/j —> JXa — > 77r"7r+7r^ is calculated as: 

Bri^' 7Xci in'aU980) ^ 7-°.fo(980) ^ 7-°-+-") = ^("0 /o) Aa/ + iV(/o)/e, ^ 

N^' ■ isr(7r" — > 77) 

where £/ = (20.5 ± 0.2)% is the efficiency for the underlying process ip' — ?> 77r°/o(980) — ^ 77r°7r+7r~. The branching 
ratio is then determined to be (2.7 ± 4.2) x 10^^, where the error is statistical only. 



TABLE II: Results of fits to the tt+tt invariant mass spectrum in the Xa region. The fit is described in the text and the 
yield error is only statistical. N(ao — )■ /o) is the number of mixing events ao(980) /o(980). S(ao — )■ /o) is the significance 
of the mixing signal ag(980) ./o(980). N(/o(980)) is the number of /o(980) events from other processes. S(/o(980)) is the 
significance of the /o(980) contribution from other processes. 



Mode 


N(a« ^ /o) 


S(aS ^ /o) 


N(/o(980)) 


S(/o(980)) 


mixing + /o(980) + 2nd poly. 


6.4±3.2(< 11.9) 


1.9cr 


0.0±8.6(< 16.7) 


< O.lcr 


mixing + 2nd poly. 


6.4±3.2(< 12.1) 


3.0a 






/o(980) + 2nd poly 






12.8±6.7(< 23.6) 


2.3a 



If we fit the TT~^Tr~ invariant mass spectrum only with the mixing signal plus a 2nd order polynomial background, 
the fit yields 6.4 ±3.2 events for the mixing signal. Comparing with the fit result with only the 2nd order polynomial, 
the change in —InL with A{d.o.f.) = 1 is 4.41, corresponding to a statistical significance of 3. Oct. The upper limit at 
the 90% C.L. is 12.1 events. 

If we assume there is no mixing and fit the tt+tt^ invariant mass spectrum only with the /o(980) contribution from 
other processes plus a 2nd order polynomial, the fit yields 12.8 ± 6.7 events for the /o(980) contribution from other 
processes. Comparing with the fit result with only the 2nd order polynomial, the change in —InL with A(d.o./.) = 1 
is 2.62, corresponding to a statistical significance of 2.3ct. The upper limit at the 90% C.L. is 23.6 events. The fit 
results are listed in Table II. The total branching ratio of ■0 — > 7Xci ~^ 71"° /o (980) — > 77r°7r+7r^ is calculated to be 
(6.0 ± 3.1) X 10"^ where the error is statistical only. 
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FIG. 4: Results of the simultaneous fit of the ■k'^tv~ mass spectra: (a) in the Xci mass window, (b) in the Xci sideband region. 
The solid curve is the result of fit described in the text. The dotted curve shows the mixing signal. The dash-dotted curve 
indicates /o(980) from other processes. The dashed curves in (a) and (b) denote the background polynomial. 

VI. DISCUSSION 



Various models for the ao(980) and the /o(980) [2-6] give different predictions for their coupling constants and 
masses; these have been measured by several experiments [19, 21-25]. From these theoretical and experimental values 
of the resonance parameters, predictions for ^af and ^/a are calculated [13, 14]. Using the parameter sets listed in 
Table III, the line shapes of the 0(5(980) -> /o(980) and /o(980) ao(980) mixing signals can be determined from MC 
simulation and the underlying /o(980) (ao(980)) shapes can be parameterized accordingly. Table III shows the fitting 
results obtained using a similar fitting procedure as described in Sects. IV and V. The fitting results are consistent 
within the statistical error. 



VII. SYSTEMATIC UNCERTAINTIES 

The systematic uncertainties on the branching ratios are summarized in Table IV. 

The systematic uncertainty associated with the tracking efficiency has been studied with control samples such as 
J/tp pTT, J/ip ppTT^TT^ , J/ip ^ K*K KsKtt. The difference of the tracking efficiencies between data and MC 
simulation is 2% per charged track. 
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TABLE III: Fitting results with various theoretical and experimental values of the resonance parameters. 



Mixing shape 


N(/o ^ al) 


S(/o ^ al) 


N(ag ^ /o) 


S(ao ^ h) 


qq [2] 


19.8±8.6(< 31.8) 


2.4a 


5.9±2.8(< 10.9) 


2.0a 


q'<f [2] 


19.4±8.5(< 31.5) 


2.4a 


6.3±3.0(< 11.5) 


2.1a 


KK [3-5] 


14.5±10.8(< 28.3) 


1.3a 


5.8±2.7(< 10.5) 


1.6a 


qqg [6] 


25.4±9.7(< 38.2) 


2.9a 


6.6±3.2(< 12.2) 


2.6a 


SND [21, 22] 


21.7±9.3(< 33.1) 


2.5a 


6.0±2.9(< 11.1) 


2.0a 


KLOE [23, 24] 


23.3±8.0(< 34.9) 


3.3a 


6.3±3.0(< 11.6) 


2.0a 


BNL [25] 


28.7±6.8(< 38.7) 


4.1a 


6.4±3.0(< 11.8) 


2.4a 


CB [19] 


27.1 ±8.4(< 37.8) 


3.7a 


6.4±3.1(< 11.8) 


2.2a 



The uncertainties due to PID of tt and K are determined from studies of control samples such as J/i/' — > /ott, 
J/'4) — > pp7r+7r~, J/ij] ^ K*K K^K^ir^ . The difference of the PID efficiency between data and MC is 2% per 
track. 

The uncertainty due to photon detection and photon conversion is 1% per photon. This is determined from studies 
of the photon detection efficiency in control samples such as J /ij} — > p'^i:^ and a study of photon conversion via 

e+e^ — > 77. 

The uncertainty due to the tt*^ selection is determined from a high purity control sample of J/?/' tt+tt^tt'^ 
decays. The tt*^ selection efficiency is obtained from the change in the tt*^ yield in the tt+tt" recoiling mass spectrum 
with or without the tt*^ selection requirement. The difference of the tt*^ reconstruction efficiency between data and 
MC simulation gives an uncertainty of 2.0% per 7r°. The uncertainty from the 77 selection is 2.0% per 77, which is 
determined in a similar way from a control sample of J /i}} — !• pprj. 

The uncertainty of the kinematic fit for J/-0 -> (/)/o(980) — > 0ao(98O) — > ^jyTr" is estimated from J/?/' — > a;?7 — s- 
TT'^^TT^Tr'^r; (r/jTr'^ 77). The efficiency is obtained from the change in the yield of uj signal by a fit to the tt^tt^tt'^ 
mass distribution with or without the requirement of the kinematic fit. The systematic uncertainty is determined to 
be 0.9%. The uncertainty of the kinematic fit for tp' 7Xci — > 77r'^ag(980) — > 77r'^/o(980) — !■ 77r'^7r+7r^ is estimated 
to be 1.7% from ip' tt+tt^ J/t/j, J/^|: -> 777. 

The branching ratios for the 77 — > 77, 7r° -> 77 and <f> ^ K^K^ decays are taken from the PDG [17]. The 
uncertainty on these branching ratios are taken as a systematic uncertainty in our measurements. 
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The total number of J /ij: events is (2.252 ± 0.028) x 10®, determined from inelusive J/ij] hadronic decays [15], and 
the total number of -if)' events is (1.06 ±0.04) x 10^, determined from inclusive ip' hadronic events [16]. The uncertainty 
on the number of J/i/' events is 1.3%, and the uncertainty on the number of ip' events is 4%. 

To estimate the systematic uncertainties due to the fit procedure, we repeat the fit with appropriate modifications 
to estimate the systematic uncertainties. The largest difference of the yield of each sources with respect to the values 
derived from the standard fit is considered as a systematic error. We change the sideband range and the order of 
polynomial to estimate the uncertainty from the background shape. A series of fits using different fitting ranges is 
performed and the largest change of the branching ratios is assigned as a systematic uncertainty. 

The total systematic uncertainties for the branching ratio measurements are obtained by adding up the contributions 
from all the systematic sources in quadrature. 

The uncertainty due to the parameterization of the mixing signal line shape and the underlying ag(980) (/o(980)) 
is kept separate and quoted as a second systematic uncertainty. The uncertainty is obtained by comparing the results 
with the parameter sets in Table III with the standard fit. We take this difference as a conservative estimate of the 
uncertainty and assign an uncertainty of 43.8% for the a[J(980) /o(980) mixing measurement and an uncertainty 
of 9.4% for the /o(980) — > ao(980) mixing measurement. For the total branching ratio measurement of J/ip 
(^00(980) <j>7]Tr'^, the uncertainty is assigned to be 38.3% and for the total branching ratio measurement of -0' 
7Xci 77r°/o(980) — > jir'^n^TT^, the uncertainty is assigned to be 9.4%. If we assume there is no mixing, the 
uncertainties of the total branching ratio measurements are assigned to be 41.4% and 42.2%, respectively. 

VIII. SUMMARY 

Based on (2.252 ± 0.028) x 10® J/ip events and (1.06 ± 0.04) x 10® ip' events, the mixing branching ratios are 
measured to be: Br{J 0./o(98O) ^ <?!)aj^(980) ^ (pr^TT^) = (3.3 ± 1.1 (stat.)±0.4 (sys.)±1.4 (para.)) x 10~^ 
Br(V'' ^ 7Xci ^ 77r0a[](980) 77r"/o(980) 77r"7r+7r-) = (2.7 ± 1.4 (stat.)±0.7 (sys.)±0.3 (para.)) x 10"^ where 
the uncertainties are statistical, systematics due to this measurement, and systeniatics due to the parameterization, 
respectively. 

The total branching ratio of J/-0 — > (^ao(980) cprjir'^ is measured to be (5.0 ±2.7 (stat.)±1.7 (sys.)±1.9 (para.)) x 
10~^ and the total branching ratio of -0' ^ IXci ^ 77r°/o(980) — > 77r°7r+7r~ is measured to be (2.7 ± 4.2 (stat.)±3.6 
(sys.)±0.3 (para.)) x 10^^. If we assume there's no mixing, the total branching ratios are measured to be (9.7 ± 2.2 
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TABLE IV: Summary of systematic errors for the branching ratio measurements and upper limits determination. 



J/V' ^ (^aJl(980) <j>rjn° 


mixing Br 


upper limit of mixing Br 


total Br 


total Br(no mixing) 


Charged tracks 


4.0% 


4.0% 


4.0% 


4.0% 


Photon detection 


4.0% 


4.0% 


4.0% 


4.0% 


PID 


4.0% 


4.0% 


4.0% 


4.0% 


77 construction 


2.0% 


2.0% 


2.0% 


2.0% 


Tv^' construction 


2.0% 


2.0% 


2.0% 


2.0% 


Kinematic fit 


0.9% 


0.9% 


0.9% 


0.9% 


Intermediate decay 


0.54% 


0.54% 


0.54% 


0.54% 


Normalization 


1.3% 


1.3% 


1.3% 


1.3% 


Background shape 


6.6% 




28.7% 


4.9% 


Fitting range 


6.2% 




14.5% 


12.9% 


Total 


11.9% 


7.7% 


33.1% 


15.8% 


V'' ->7Xci ^7'r'Vo(980) 


mixing Br 


upper limit of mixing Br 


total Br 


total Br (no mixing) 


Charged tracks 


4.0% 


4.0% 


4.0% 


4.0% 


Photon detection 


3.0% 


3.0% 


3.0% 


3.0% 


PID 


4.0% 


4.0% 


4.0% 


4.0% 


■k'^ construction 


2.0% 


2.0% 


2.0% 


2.0% 


Kinematic fit 


1.7% 


1.7% 


1.7% 


1.7% 


Intermediate decay 


0.034% 


0.034% 


0.034% 


0.034% 


Normalization 


4.0% 


4.0% 


4.0% 


4.0% 


Background shape 


23.4% 




128.1% 


48.4% 


Fitting range 


6.3% 




32.8% 


17.2% 


Total 


25.5% 


8.0% 


132.5% 


52.0% 



(stat.)±1.5 (sys.)±4.0 (para.)) x 10"^ and (6.0 ±3.1 (stat.)±3.1 (sys.)±2.5 (para.)) x lO^'^, respectively. 

When determining the upper limit of the number of J/'i/j — > 0/o(98O) — !> 0ag(98O) — > (f>7]Tr'^ events, the uneertainties 
due to the fit range, background shape, the parameterization of mixing signal line shape and the underlying ao(980) 
are considered. Using the Bayesian method, different upper limits at the 90% C.L. are determined by varying the fit 
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range, the background shape and the parameterization of the mixing signal hne shape and the underlying ag(980) in 
Table III. The upper limit for the mixing signal is taken to be the largest of them: N^^^' =39.7 events. A conservative 
estimate of the upper limit of the branching ratio is determined by lowering the efficiency by one standard deviation. 

The upper limit on the branching ratio at the 90% C.L. is calculated as : 

Br{J/i} 0/o(98O) ^ (/)aj](980) ^ 0r;7r°) 

^U.L. 

<^ L2 = 5 4 X 10^^ 

{efa-<jf,^-)-Nj/^-Br{c^^K+K-)-Br(7^^^-f)-Br(TT'' ^-f-f) 

Similarly, considering the uncertainties due to fit range, background shape, the parameterization of mixing signal 
line shape and the underlying /o(980) , the upper limit number of the mixing signal ip' — > 7Xci ^ 77r*'aQ(980) 
77r°/o(980) 77r"7r+7r~ is determined to be: N^j^ = 13.0 events. 

The upper limit on the branching ratio at the 90% C.L. is calculated as : 

Br{il)' 7Xci ^ 77r°a[J(980) ^ 77r"/o(980) ^ 77r"7r+7r~) 

< Tj-r- — r-?^ r = 6.0 X lO"'^. 

(-/-^^/■)-^^^'-Sr(7r"^77) 



The mixing intensity for the /o(980) a§(980) transition is calculated to be: 

_ BrjJ/iP ^ 0/o(98O) ^ 0a°(98O) ^ 0??7r") 
" Br(J/V' ^ 0/o(98O) ^ 07r7r)[2o] 

= (0.60 ± 0.20{stat.) ± 0.12(sys.) ± 0.26(para.)%, 

where the uncertainties are statistical, systematics due to this measurement, and the parameterization, respectively. 
The uncertainty from Br{J — > ^/o(980) — > ^tttt) is included as a part of the systematic error. The upper limit of 
the mixing intensity at the 90% C.L. is 1.1%. 

The mixing intensity ^a/ for the ao(980) — > /o(980) transition is calculated to be: 

_ Br{i]y 7Xci ^ 77r0a0(980) -> 77r°/o(980) ^ -^-k^tt+'k-) 
~ Brill;' 7Xci -> 77r"a[](980) ^ 77r07rO?7)["l 

= (0.31 ± 0.16(stat.) ± Q.lA{sys.) ± 0.03(para.))%. 

The uncertainty from Br{ip' — ^ 7Xci 77r°aQ(980) -> 77r*'7r°77) is included as a part of the systematic error. The 
upper limit of the mixing intensity ^af at the 90% C.L. is 1.0%. 
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The calculated mixing intensities [14] with the resonance parameters from various models [2-6] and experimental 
measurements [19, 21-25] are compared with our results in Fig. 5. This result will be very useful in pinning down the 
resonance parameters of ao(980) and /o(980). 




10" lO-" 

mixing intensity ^ 



fa 



FIG. 5: Mixing intensities of /o(980) aS(980) and a8(980) /o(980).Tlie dots are predictions for the mixing intensities 
with theoretical and experimental values of the parameters [14]. The shaded region is our results with statistical errors and 
systematics due to this measurement and the parameterization. The solid lines mark the upper limits. 
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